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We study the non-equilibrium dynamics of a one-dimensional system of hard core bosons in the 
presence of an onsite potential (with an alternating sign between the odd and even sites) which 
shows a quantum phase transition from the superfluid (SF) phase to the Mott Insulator (MI) phase. 
The ground state quantum fidelity shows a sharp dip at the QCP while the fidelity susceptibility 
shows a divergence right there with its scaling given in terms of the correlation length exponent of 
the quantum phase transition. We then study the evolution of this bosonic system (with a qubit 
globally coupled to it) following a quench in which the magnitude of the alternating potential is 
changed from zero (the SF phase) to a non-zero value (the MI phase) according to a half Rosen Zener 
scheme. The loss of coherence of the qubit (initially in a pure state) after the quench is investigated 
by calculating the time dependence of the decoherence factor. This result is compared with that 
of the sudden quench limit of the half Rosen-Zener scheme where an exact analytical form of the 
decoherence factor can be derived. The local von Neumann entropy density is calculated in the final 
MI phase and is found to be less than the equilibrium value (log 2) due to the defects generated in 
the final state as a result of the quenching that starts from the QCP of the system. We also briefly 
dwell on the full Rosen-Zener quenching scheme in which the system is finally brought back to the 
SF phase through the intermediate MI phase and calculate the reduction in the supercurrent and 
the non-zero value of the residual local entropy density in the final state. 



I. INTRODUCTION 

Recent advancements in experiments on ultracold 
atoms trapped in optical lattices have facilitated the re- 
alization of ultracold vapors of bosonic atoms, and hence 
have opened up new directions towards the experimen- 
tal studies of low dimensional bosonic systems^ ^. For 
example, following the pioneering experiments indicat- 
ing a superfluid (SF) to a Mott insulator (MI) transition 
in optical lattices in three-dimensiorF\and also in one 
dimension^) and the corresponding study on the non- 
equilibrium dynamics^, there is an upsurge in the studies 
of quantum phase transitions (QPTs)^ and dynamics 
of trapped atoms in optical lattices. More interestingly, 
two dimensional optical lattices have made the quasi 
one dimensional regime experimentally accessible^ by 
keeping the transverse potentials much higher than the 
longitudinal potential. By appropriately tuning the lon- 
gitudinal potential, different limits of the bosonic Hub- 
bard model has been realized. One of such limits happens 
to be the hard-core boson (HCB) limit (or the Tonks- 
Girardeu^^ limit), where two bosons can not occupy 
the same s ite; th is limit has also been achieved in an op- 
tical lattice^^n^. These experiments have paved the way 
for a plethora of the oretical studies in low-dimensional 
bosonic systems^l^^ e speci ally from the viewpoint of 
the SF t o MI t ransitiod-^^^ and related non-equilibrium 
dynamic d^^ l ^-^ i The HCB systems have turned out to be 
very advantageous in this contextP^H^. 

In parallel, there have been numerous studies which 
attempt to bridge a connection between QPTJ^H^I 
and quantum information theoretic measures like 
concurrence^^ entanglement entrop}/^^ quantum 
fidelit}^^, quantum discord^^, etc.. These measures 
enable us to detect a QCP and they also show distinctive 



scaling relations close to it characterized by some of the 
associated critical exponents. Similarly, the decoherence 
(or loss of phase information)^ of a qubit coupled to a 
quantum critical system is also being investigated^^. 

The scaling of the density of defects (or heat) pro- 
duced following a sIoa^^^^ or rapid quenching^^ across 
(or starting from) a QCP has also attracted attention of 
the scientists. Defects generated in the final state of the 
quantum system due to the quenching through a QCP in 
turn lead to non-zero quantu m corr elations (for example, 
non-zero local entropy densit}0^Ell, concurrence^, quan- 
tum discord^^, etc.) in the final state which are otherwise 
absent in the defect free final state. These information 
theoretic measures have also been found to satisfy scaling 
relations identical to that of the defect density in some 
cases. For recent reviews, see [49-51 . 

In this paper, we study the dynamics of a one- 
dimensional lattice of HCBs in which Bosons are sub- 
jected to an onsite potential. The model has a SF long- 
range order which persists up to a threshold value of the 
onsite potential at which there is a QPT from the SF 
to the MI phase. We are however interested in the case 
where the onsite potential is site-dependent (rather, al- 
ternates in sign on the even and odd sites); under this 
condition the SF long-range order is destroyed as soon as 
the potential is switched on. We note that this model has 
been studied under a half Rosen-Zener (HRZ) quenching 
schem#^^ in which the magnitude of the alternating on- 
site potential is quenched from zero to a non-zero value 
and the residual supercurrent in the MI phase has been 
estimatecP^l. 

The motivation of this work is the following: although 
there has been a series of studies of quantum critical 
dynamics which involve Landau-Zener tunneling^^ (for 
many examples, see j49H5T] ), the Rosen-Zener (RZ) tun- 
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neling (for which the non-adiabatic excitation probabil- 
ity can also be exactly calculated) has received rela- 
tively less attention. We use the integrability of the 
one-dimensional HCB system in an alternating poten- 
tial along with the exact analytical results for the HRZ 
quenching to investigate the decoherence of a qubit con- 
nected to the HCB system following a HRZ quenching of 
the magnitude of the onsite potential. We also calculate 
the generation of local entropy in the HCB system in its 
final MI state following the quench. Given the current 
interest in QPTs, dynamics and quantum information as 
discussed above, these results are expected to be useful 
both from experimental and theoretical viewpoints. 

The paper is organized in the following way: in Sec. 
H, we describe the QPT in the HCB chain in an alter- 
nating potential by analyzing the energy spectrum of the 
Hamiltonian; any non-zero value of the alternating po- 
tential leads to an energy gap in an otherwise gapless 
spectrum so that the system is in the MI phase. In Sec. 
HI, we show how this QPT can be detected and char- 
acterized by investigating the ground state fidelity and 
fidelity susceptibility. 

The dynamics of the HCB chain is studied in Sec. IV. 
Here a qubit (or a central spin- 1/2) is globally coupled to 
the HCB chain. Our focus is limited to the case when the 
coupling between the qubit and the HCB chain, which in 
fact plays the role of an environment to which the qubit 
is coupled, is very weak. We study the decoherence of 
the qubit by calculating the decoherence factor in the fi- 
nal state when the onsite potential is changed from zero 
(the SF phase) to a finite value (the MI phase) following 
a HRZ quenching scheme. An exact expression of the 
decoherence factor is derived analytically in the sudden 
quench limit where the alternating potential is instan- 
taneously switched on and the results are compared to 
those of the previous case. 

Finally, in Sec. V, we investigate the single site (lo- 
cal) von Neumann entropy density in the final MI phase 
following the HRZ quenching. We note that the local en- 
tropy density is zero in the SF phase and is equal to log 2 
in the MI phase because of its bipartite structure. We, 
however, find that the value of this entropy in the final 
MI phase reached after the quenching is less than log 2 
by an amount which depends on the parameters of the 
HRZ quenching. This deviation is due to the fact that 
the system is quenched out of the SF phase (which is also 
a gapless QCP) at a finite rate which leads to the defects 
that in the final MI phase leading to a surviving super- 
current and reduced local entropy density. In Sec. VI, we 
study the HCB chain under the full Rosen Zener (FRZ) 
quenching scheme in which the system is finally brought 
back to the SF phase through the intermediate MI phase 
and the surviving supercurrent and the residual local en- 
tropy density are calculated. 



II. THE MODEL 

We consider the Lattice-Tonks-Girardeu gas (hard- 
core) limit of the one-dimensional Bosonic Hubbard 
model^^ given by the Hamiltonian 

H = -w Y.(h\hi+i + h.c) + V Y,{-^y^% (1) 

I I 

where w is the hopping amplitude, V is the onsite po- 
tential; hi and h\ are the bosonic annihilation and cre- 
ation operators at the l^^ site of the lattice, respectively. 
These bosonic operators satisfy the canonical commuta- 
tion relation \b\^hm] = ^im] additionally, the hard core 
condition demands, (6^)^ = = {b])'^. The Hamiltonian 
([T]) undergoes a QPT from the gapless SF phase to the 
gapped MI phase for any non-zero value of the alternat- 
ing potential V as shown below. 

This Hamiltonian can be exactly solved using Jordan- 
Wigner (JW) transformations^^ given by 



4 = 



Yl exp {al^ar, 



.m<l 



(2) 



where a] and ai are the JW fermionic operators satisfy- 
ing the fermion anti-commutation relations {aj^am} = 
^im^{(ihCim} = 0. Using JW transformation followed 
by the Fourier transformation, the energy spectrum of 
Hamiltonian ([T]) can be exactly obtained. In terms 
of JW fermions, the Hamiltonian can be re- written as 
H = Ho ^Hdj where. 



Ho = - ^ 2w cos k{alak - ^^I+t^^/c+tt), 

|/c|<7r/2 
|fe|<7r/2 



(3) 



Evidently, the mode with wave vector k couples to the 
(A: + tt)— mode, one can rewrite the Hamiltonian in the 
reduced 2x2 form. 



n = 



|/c|<7r/2 



with 



2w cos k 
-V 



-V 
-2w cos k 



and the energy spectrum (see Fig.Q) is given by 



Ek = V'^w^ cos2fc + T/2_ 



(4) 



(5) 



(6) 



We note that the spectrum (|6| is gapped even for an 
infinitesimal alternating potential implying that the sys- 
tem is in the MI phase for any V 7^ 0. On the other 
hand, for V = 0, the spectrum is gapless for the critical 
mode k = 7r/2, and the spin chain is in the SF phase. 
From the spectrum, we find that the QPT at F = is 
characterized by the correlation length exponents v = 1 
and the dynamical exponent z = 1. 
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FIG. 1: (Color online) The energy spectrum ([6| for the model 
|T]) in the reduced Brillouin zone. In the SF phase = 0), 
the spectrum is gapless at /c = 7r/2 (red dashed). A non-zero 
^generates a gap in the excitation spectrum at the critical 
modes are at /c = ±7r/2. 



where we have assumed a small system size (L) and also 
S ^ limit, which allow us to truncate the above series 
at the S'^ order; in the present problem spatial dimension- 
ality d = 1. The quantity xf = -(2/L^n(F)/j2|^^^^ 
called the fidelity susceptibility densit}^!^'^, is a measure 
of the rate of the change of the ground state wave function 
when the parameter V is changed infinitesimally. Usu- 
ally quantum fidelity shows a sharp dip at a QCP where 
Xf diverges with the system size; the universal scaling 
of Xf is given in terms of some of the critical exponents 
associated with the QPT. 

To calculate F and xf in the vicinity of the QPT of 
Hamiltonian ([T]), we use the reduced two-level Hamil- 
tonian ([5|. One can use Bogoliubov transformation to 
obtain the ground state wave function for a particular 
momemtum mode and a given potential V in the form 



\Mk. V)) = cos{Ok{V))\k) + sm{Ok{V))\k + n) (8) 



III. FIDELITY AND FIDELITY 
SUSCEPTIBILITY 

One of the most widely used quantum information 
theoretic measure for detecting and characterizing quan- 
tum phase transitions is ground state quantum fidelit}/^^ 
which is the magnitude of the overlap of the two ground 
states of a quantum many body system belonging to dif- 
ferent values of a parameter of the Hamiltonian. Re- 
ferring to the Hamiltonian ([T]), we can define quantum 
fidelity F{V^ V -\- S) between two ground states with the 
alternating potentials V and U + (5, respectively, given by 



{Mv)\Mv + s)) 



XF 



(7) 



where tan(2^/e(U)) = —V/{2wcosk). An exact expres- 
sion of quantum fidelity can be then obtained using 
Eqs. ([7| and ([8|: 



F = Y[\cos{0k{V)-0k{V^5)\ 



exp 



2n 



dk\cos{ek{v)-Ok{v + s)\ 



TT — TV/L 



■ (9) 



Expanding around the critical mode k = 7r/2, one ar- 
rives at the simplified form 



F = exp 



32aV^ 



|tan "^(a)— tan ^[a{L — l) 



-1)2/ 



a = 



2W'K 

Vl 



(10) 



As shown in Fig.([2|, the fidelity shows a dip and the 
susceptibility shows a peak at the QCP, V = 0. We also 
find that xf scales as L near U = (SF phase) and 
as deep inside MI phase (see Fig.pl) . This is in 
congruence with the generic scalin 

near the QCP {L <C U"^), and xf T^""^"^ away from 
the QCP (L > y-^), with u = d = 1. 



IV. DECOHERENCE AND HRZ QUENCHING 

A. HRZ quenching for the HCB chain with a qubit 
coupled to it 

In this section, we shall explore the decoherence of a 
qubit coupled to the environment, chosen to be the HCB 
chain ([T])), which is driven following a quenching scheme. 
We assume a global coupling between the qubit and all 
the bosons of the model ([T]) with the coupling Hamilto- 
nian given by 



(11) 
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FIG. 2: (Color online) Fidelity shows a clear dip at V = 0. 
The inset shows that a peak occurs in the fidelity susceptibil- 
ity at the critical point. 
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FIG. 3: (Color online) Numercially obtained scaling oi xf- 
away from the QCP, Xf while in the vicinity QCP, 

Xf ^ L (see inset). These scaling relations are in agreement 
with the theoretical prediction. 



channels, denoted by l-L^ and , respectively, are given 
by 

±_f2wcosk -{V±S)\ . . 

- y_(^v±S) -2wcoskJ' ^^^^ 

We shall denote the corresponding time-evolved states of 
the environmental Hamiltonian corresponding to these 
two branches as |V^+(t)) and \ilj~{t)), respectively. 

In this paper, we shall employ the HRZ quenching 
scheme in which the alternating potential is changed from 
zero to a finite value Vb, (see Fig. Q) in a non-linear 
fashion given h^^^^ 

This implies that the system is quenched from the SF 
phase {t — oo) to the MI phase; our aim is to analyze 
the time evolution of the environment and the resulting 
decoherence of the qubit coupled to it in the MI phase 
for t > 0. One should also note that the spectrum (see 
Eq. (|6|) of the final Hamiltonian depends on Vb but not 
on r. This implies that Vb and r are not on the same 
footing for the HRZ quench as will be reflected in results 
obtained below; this is in contrast to the FRZ quench 
discussed in Sec. VI. 



VOM _ 

/ / Ti > T2 > T3 



where cr| represents the qubit, bjbi is the number density 
of the environmental HCB chain at site 1; S is the coupling 
parameter between system and environment. (The form 
of the coupling Hamiltonian (11) can be interpreted in 



the following way: the HCB chain can be recast to a 
transverse XY spin chain in a transverse field in the z- 
direction; the z component of the spin at the site / is 
coupled to the z-component of the central qubit.) In 
subsequent sections. We shall work in the limit of a weak 
coupling between the central qubit and HCB system (i.e.. 

Due to the coupling to the central qubit, the time evo- 
lution of the environmental bosonic chain is split into two 
channels, corresponding to the 1 1) (= +1) and 1 1) (= —1) 
state of the the qubit. Using Eq. ([5|, we find that the 
reduced Hamiltonians of the HCB system for these two 



FIG. 4: (Color online) The HRZ quenching scheme for V{t). 
We get the sudden quench limit for r ^ 0. 

One can show that in the limit (5^0, the off-diagonal 



terms of the Hamiltonian( 12 ) can be written as 



V{t)±S = V^{t) 



Tit 

(Vb ± ^)sech^. 



(14) 



where = r ±a^r, with a being a number of the order 
of unity. When compared with the RZ form Eq. (13), 



this approximation implies the following: in the limit of 
a very weak coupling between the qubit and the environ- 
ment, the evolution of the two channels can be viewed as 
two independent HRZ quenches with final potentials and 
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quenching parameters [(Vq + 6)^ r+] and [(Vq ~ ^~]^ 
respectively. 

In order to calculate the time evolution of 
and \ip~{t)) at a given instant t, let us consider a 
generic state for a given momentum mode: \i^^{t)) = 
s^{t)\k) -\- p^{t)\k + tt). Using Schrodinger equations 
id\ip^ (t)) / dt = y^^\i^^{t)) ^ it can be shown that time 
evolution of the probability amplitudes s^{t) and p^{t) 
are dictated by the equations, 



^{t)2w cos{k)^p^{t)V^{t), 
-p^{t)2wcos{k)^s^{t)V^{t). 



(15) 



Using transformations S^{t) = ex.p {2i wt cos k)s^{t), 
P^{t) = ex.p{—2iwt cos k)p'^ (t) ^ and utilizing the form 
of V^{t) from Eq. ([l4|, we get 



5± 



^cosh(f^) 
Aiw cos{k) - 



(16) 



-tanh 



TTt 



which can be reduced to a hypergeometric form with the 
initial conditions, (X))| = 1 and |P^(— oo)| = 0. 

Expanding near the critical mode {k = 7r/2), one even- 
tually finds the solution at t = of the form: 



P^(0) 



s±(0) = cos 



(17) 



Exploiting the continuity condition of the wave func- 
tion at t = 0, let us write the generic wave function for 
t > in the form 



it)) 



|5^(0)) 



c^e 



where |^^) and |e^) are the ground state and excited 
state wave functions (with energies -E^ and E^) with 
probability amplitudes and c^, respectively. Express- 
ing Eq. ([T8| in terms of momentum modes \k) and \k-\-7r) 
and using Bogoliubov transformation, we get 



-p^m^)\k) 

-p±(0)Af)|fc- 



(19) 



where A^{t) 



cos{E^t) 



zcos{20^)sm{E^t),B^{t) = 
ism{20^)sm{E^t) and E^ = y/4w^cos^ k^ {Vo±S)^. 
In the next section, we shall use Eq. (IT9|) and analyze 



the behavior of decoherence factor of the qubit coupled 
to the HCB chain. 

B. Decoherence of the qubit 



assume that the qubit is initially in a pure state at 
t ^ —oo. The off-diagonal terms of the reduced den- 
sity matrix for t > incorporate the decoherence factor 
D{t) = \{ilj~^{t)\ip~ {t))\'^ ^ which measures the decoher- 
ence of the qubit. A non-zero value (less than unity) of 
D{t) implies that the qubit is in a mixed state and initial 
phase coherence is lost; D{t) = denotes perfect mix- 
ing. Considering the two-level structure of the reduced 
Hamiltonian of the environmental HCB chain (see (12)) 
we get. 



Dit) = n \Dkit)f; D,it) = mmkim, (20) 

k 

which can be put in the form 

D{t) = \{^+m-m' 



tt/2 



exp ( 1^ / \n{\Dk{t){')dk 

(21) 

To evaluate -D(t), we now use Eq. (19) and work in the 
limit (5 0, when one can approximate ff^ (defined after 
Eq. (8) with F ^- Vo ± (5) as 



^t 



S — — 
05 



(22) 



where 



dot 
05 



5=0 



dS 



6=0 



-2w cos k 



6=0 



Aw'^ cos^ k^V^;' 



(23) 



One can obtain using Eq. (23) 



(0)), (18) {i,t{mkit))=^os{Ett-E^t)cos{2^)+ 



do 

sm{E^t - E't) sin(27) ( sin 2l9fc + 2 cos 2l9/e 5^ 



\ sm{E^t - E't) cos(20) ( cos 20^ + 2 sin 20^ S 



86 



6=0 



6=0 

(24) 



where 7 = St{1 + Voa)/2 and ^ = Vbr/2. 

The maximum contribution to the Eq. ( [24| ) comes from 
the modes close to the critical mode k = ttJ2. We assume 
small r (r ^ 1) so that Sr ^ and use the fact that 
ln(l — x) ^ —x; the decoherence factor in the early time 
limit gets reduced to the form 



2 i2 



4^2^/2 ^ y2 



(l^Ar^w^k'^), (25) 



To study the decoherence of the qubit coupled to the 
HCB chain following the HRZ quench (t > 0), one in- 
vestigates the reduced density matrix of the qubit. We 



where k' = 7t/2 - k. Using Eq. ^ and Eq. ^ and 
retaining terms up to the leading orders in S and r, we 
get 



6 



InD(t) 



L6H^ r 2Vo 
27r \ w 



tan 



Vo'r 



(26) 




time log{Vb) 



FIG. 5: (Color online) (a) Gaussian fall of D(t) following a 
HRZ quenchis shown with time. The different lines corre- 
spond to different values of Vq. Increasing Vb makes the de- 
coherence fall faster. The range of Vb covered in the plot is 
from 1 (blue) to 50 (black). For Vb > ^ttw, D{t) is approxi- 
mately independent of Vb which is refelcted by the bunching 
up of the curves for different Vb for higher values, (b) The 
similar nature is observed for D(t) in the sudden quench case 
(r = 0). 



Eq. ( 26 ) is plotted in Fig. ([5| which shows a gaussian 
fall in time of the decoherence factor in the early time 
limit; this Gaussian fall is the expected behavior in the 
vicinity of a QCF^^. To analyze the scaling of the de- 
coherence factor with Vb near the QCP (Vo = 0), we 
define a quantity A(Fo) = — log^(^); Fig. Q shows that 
^4 ^ Vq. On the other hand, when Vq exceed a threshold 
value {=27rw)^ one can expand the tan~^ term in Eq. (26) 



and show that D{t) is approximately independent of Vq. 
Moreover, Fig. ([5| also shows that logl^(t) depends very 
weakly on r for r <C 1, which is further illustrated in 
Fig. (|6| where we show that curves for different values of 
r fall on top of each other. 



C. Decoherence in the sudden quench limit 

The time evolution of decoherence factor can be de- 
rived exactly from the overlap of the initial wave function 
and the time-evolved version of the final wave function 
in the sudden quench (r = 0) limit of the HRZ scheme 
as shown below. In this limit, the onsite the potential 
is y = for t < and it is abruptly changed to Vq -\- S 
{or Vq — S depending on the initial state of the qubit) at 
t = 0. The ground state of the HCB system at t = is 



FIG. 6: (Color online) Numerically obtained scaling: log D ^ 
Vq. It can be seen that the curves for different values of r 
coincide with each other highlighting the weak dependence of 
D{t) on r. 



|?/^(0)) = which can be written with the Bogoliubov 
parameters as 



|fc)=cos^±|ff±(0))-sin^±|e±(0)) 



(27) 



Since the ground state and the excited state evolve in 
time with the corresponding eigenenergies, the wave func- 
tion for t > is simply given by 



(28) 

Expressing |^^(0)) and |e^(0)) in terms of \k) and \k-\-7r) 
using the Bogoliubov transformations, we find 

\^kit))sQ = |fc)[cos(S±t) +isin(E±t)cos(2e±)] 
+ \k + TT)[ism{E^t)sm{20f)], 

so that using Eq. p9) ), we find 

{i'UtMk(t))sQ = cos{E+t - E^t)- 



sm{Elt - Ef, t) jcos 26'fc + 25 sin 26'fc -^1^=0 



so that we find 



log D{t) 



27r 



(30) 



(31) 



Comparing Eq. (26) with Eq. (31), we conclude that the 



decoherence factor for the HRZ quench in the small r 
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limit has additional correction terms of the order of Vqt'^. 
It can also be shown that for Vo/2w > ({^ 0.75), the 
decay constant dictating the Gaussian decay of the deco- 
herence factor in the early time limit is greater than the 
sudden quench case in comparison to the HRZ case with 
small r. 



V. VON NEUMANN ENTROPY AND 
DIAGONAL ENTROPY OF THE HCB CHAIN 
FOLLOWING THE HRZ QUENCH 

In this section, we shall calculate the von Neumann 
entropy given by — Tr plog(p) where p is the density ma- 
trix of the final state of the system. Ideally in the MI 
phase, the local von Neumann entropy density s = log 2. 
(The MI phase is in a pure state and hence the global 
entropy is zero. However, the (single site) local entropy 
obtained by integrating over the momentum modes is 
non-zero because of the bipartite structure of the MI 
phase. Interpreting in terms of the spin variables, when 
observed locally upon "coarse-graining" in momentunP^ 
both the spin states appear with an equal probability 
(= 1/2) which makes the entropy density log 2 ) 

In the present context, however, the MI phase is 
reached through a non-equilibrium variation of the al- 
ternating on-site potential starting from the SF phase 
at a finite rate and hence the entropy density in the MI 
phase gets reduced. To calculate it, we decompose the 
density matrix in a direct product form, p = ^YlkPk^ 
where pk is the reduced density matrix for the k-th mode. 
Consequently, the entropy density turns to be 



7t/2 



7t/2 



dk Tr pklog(pk) 



(32) 



To calculate pk following the HRZ, we use Eq.( 19 ); taking 




FIG. 7: Variation of s with Vb in the MI phase (with w — 1). 
The inset justifies that s scales linearly with Vb for small Vb. 




the long time average of the terms of pk, we find 



FIG. 8: Variation of s with r in the MI phase with w = 1. 
The inset shows that logs scales lineally with logr with a 
slope =2. 



Pk 



/ (l+cos^ 26>fc cos(yor)) cos 26>fc sin 26>fc cos(Vor) ' 
2 2 

> cos26>fc sin26>fc cos(Vor) (l-cos^ 26>fc cos(Vor)) 
\ 2 2 



■ (33) 



Diagonalizing the density matrix, the entropy density can 
be expressed in terms of the eigenvalues X±^k as 



1 f^/^ 

7r/2 



dk A+,/e log(A+,/e) + A_,/e log(A_,fe), (34) 



whereA±^/c = ^ (1 ± cos 20^ cos{Vor)). The von Neumann 
entropy density s increases linearly with Vq for Vq < 2w^ 
and saturates to the maximum value of log 2 for higher 
values of Vq (see Fig. On the other hand, s is found 
to scale quadratically with r (see Fig. ([8|). As mentioned 
already that for a HRZ quenching, the parameters Vb and 
T are not on an identical footing which is also refiected 
in the scaling of s. 




time 



FIG. 9: The diagonal entropy density Sd is plotted against 
time shows an oscillatory interference pattern similar to that 
found in the surviving supercurrent in the MI phase after a 
similar as reported in^'^. 
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One can also calculate the diagonal entropy^ defined 
as Sd{t) = ^nPnn{t) log pnn{t) whcrc pnn{t) are the di- 
agonal terms of the density matrix obtained from Eq. ( 19 ) 
(without any time averaging). The diagonal entropy 
Sd{t) shows an oscillatory behavior (see Fig. (I9|) similar 



to the supercurrent in the MI phase following a similar 
quenclP^. The scaling of the diagonal entropy Sd with Vb 
and r is same as compared to the scaling of von Neumann 
entropy density s in both the region of Vq. 



VI. CURRENT AND VON NEUMANN 
ENTROPY STUDIES AFTER A FRZ QUENCH 
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In this section, we shall estimate the supercurrent and 
von Neumann entropy following a full RZ quench of the 
potential given by the following form: 



V{t) = Vo sech 



TTt 



-oo < t < +00; 



(35) 



the system is initially {t — oo) in the SF phase and 
finally brought back to the SF phase (as t ^ oo) through 
the intermediate MI phase. We study the time evolution 
of the system after the quenching process gets over (i.e., 
in the final SF phase). In the SF phase the reduced 
Hamiltonian is diagonal in the basis \k) and \k-\-7r) (with 
\k) (|/c + 7r)) being the ground state (excited state)). The 
wave-function of the HCB system immediately after the 
FRZ quench (which we set as t = 0) can be written as a 
linear combination of these basis states. 



mt = 0)) = V^-Pk\k) ^^/P~k\k^ tt), (36) 
where Pk is the RZ non-adiabatic transition formulsP^ 
2 f^or^ 



Pk sin 



sech [2TWCosk]. 



(37) 



The time-evolved wave-fuction at some later time t can 
readily be written as 



FIG. 10: (color online) The figure shows the variation super- 
current j against r. The inset shows the plot of log(j) versus 
log(T), which is a straight line with slope 2 



which is identical to the u dependence of the supercur- 
rent in the initial SF phase. Secondly, the supercurrent 
becomes independent of time after the full Rosen-Zener 
quench. This is due to the fact that at the final time 
HCB system reaches its eigen states. Thirdly, because of 
the passage through the MI phase starting from a QCP, 
the current in the final state is reduced from its initial 
value (at t — oo), (2k; sin i^/)7r, by a factor Vqt'^. It 
is also to be noted that in the the correction term of the 
super current is a function of the combination Vqt im- 
plying that Vb and r are on the same footing for the FRZ 
quenching. This result is numerically verified as shown 
in Fig. (pB. 



In a similar spirit, one can calculate the residual von 
Neumann entropy density in the final SF phase. The 
rapidly oscillating off-diagonal terms of the reduced den- 
sity matrix constructed from the wave function given in 
Eq. ( [38| , vanish over long time averagin ^^^ * ^^* so that 
the decohered reduced density matrix has a diagonal 
form. Calculating the local entropy density using this 
decohered reduced denstity matrix, one can show that 



IV^(t)) = ^/l - Pke~'^''^\k) + ^/Pke'^^^\k + tt), (38) ^ ^ ^o^ (see Fig.([ll|). One interesting point should be 



where = —2w cos k in the SF phase. In order to 
calculate supercurrent one has to apply a boost to the 
Hamiltonian which takes the form —w^i{e~'^^bjbi-^i + 
h.c)^^; consequently, the momentum value gets shifted 
from k to k -\- v. Expressing the super-current operator. 



m 



h.c), in the Fourier space, one 



can find its expectation value with respect to the boosted 
counterpart of the state ( 38 ) 



m = - 

TT 



7r/2 



7r/2 

2w sin V 



dk sin(/c + i/)(l-2P/e), 



2 {V^tY 



(39) 



The above result leads the following interesting ob- 
servations: in the limit of small sin(i^) ~ j ~ v 



highlighted here: the quenching through the MI phase 
generates defects in the SF phase which result in a re- 
duction of the supercurrent and a non-zero value of s 
both scaling as V^r'^. 



VII. CONCLUDING COMMENTS 

In this paper, we have studied the QPT and dynamics 
of a one dimensional HCB system in the presence of an 
onsite potential which alternates in sign from site to site. 
We have shown that the ground state quantum fidelity 
shows a sharp dip at the QCP (F = 0) indicating that 
the system is in the MI phase for any non-zero value of 
V . At the same the fidelity susceptibility is also found 
to diverge with the system size in a power-law fashion 
dictated by the critical exponent v (which is unity for in 
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FIG. 11: (color online) Variation of s with Vqt in the SF phase. 
The inset shows that s ^ Vqt'^ for small Vb and r. 

the present case). 

Subsequently, we have analyzed the scaling of the deco- 
herence factor D{t) of a central qubit which it is globally 
connected to the HCB system that is driven from the ini- 
tial SF phase to the MI phase following the HRZ quench- 
ing scheme. In the limit of a weak coupling between the 
qubit and the environmental HCB system, a threshold 
value of the magnitude of the alternating potential given 
by Vq = 27rw, is found to exist. Interestingly, the deco- 
herence factor grows linearly with Vb when Vq < 2ttw^ 
whereas for Vq > 27tw^ it turns to be independent of Vq. 
On the other hand, D{t) is found to depend very weekly 
on the quenching parameter r. This is due to the fact 
that the energy spectrum of the Hamiltonian in the MI 
phase reached through the HRZ quenching scheme de- 
pends only on Vb and not on r. In the sudden quench 
limit (r 0) an exact expression of D(t) is obtained. In 
the case of a finite but small r, the decoherence factor has 
been found to contain additional correction terms (scal- 



ing as (Vbr)^) to the same obtained in sudden quench 
limit. The more interesting observation is that we find 
the the existence an upper limit of Vo(= 1.5k;), above 
which there is a faster early time decay of the decoher- 
ence factor of the HRZ case for small r in comparison to 
the sudden quenching case. Therefore, above this upper 
limit there is a less mixing effect occurring for the qubit 
in the sudden quench case as compared to that of the 
HRZ case. 

We have also studied the local von Neumann entropy 
density and diagonal entropy of the HCB chain in MI 
phase following the HRZ quench. The von Neumann en- 
tropy density s scales linearly with Vb for small values 
of Vb (i.e., Vb < 2w) while it becomes independent of Vb 
for higher values of Vo- On the other hand, s is found to 
scale quadratically with r throughout. Interestingly, the 
von Neumann entropy density is found to be less than its 
expected value of log 2 in the MI phase. This is a conse- 
quence of the fact the system is quenched to the MI phase 
from the SF phase with 5 = (which is also the QCP) 
at a finite rate which leads to defects in the MI phase 
resulting in a surviving supercurrent^^ and reduced local 
entropy. 

Finally, we have calculated supercurrent and von Neu- 
mann entropy after a FRZ quench when the HCB chain 
is again brought back to the SF phase; interestingly the 
reduction in the supercurrent and s both scale identically 
as (Vot)'^ in the SF phase emphasizing their close connec- 
tion. It should also be reiterated that following the HRZ 
quenching there are a surviving supercurrent as well as 
a reduction in s in the MI phase. On the other hand, 
for the FRZ scheme it the other way round; one finds a 
reduction in supercurrent and a surviving s in the final 
SF phase. 
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